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Abstract We find some new exact cosmological solutions 
for the covariant scalar-tensor-vector gravity theory, the so- 
called Modified Gravity (MOG). The exact solution of the 
vacuum field equations has been derived. Also, for non vac¬ 
uum cases we have found some exact solutions with the aid 
of the Noether symmetry approach. More specifically, the 
symmetry vector and also the Noether conserved quantity 
associated to the point-like Lagrangian of the theory have 
been found. Also we find the exact form of the generic vec¬ 
tor field potential of this theory by considering the behavior 
of the relevant point-like Lagrangian under the infinitesimal 
generator of the Noether symmetry. Finally, we discuss the 
cosmological implications of the solutions. 

1 Introduction 

We study a Scalar-Tensor-Vector theory which has been in¬ 
troduced as an alternative to particle dark matter [1]. This 
theory is referred as MOG (Modified Gravity) in the current 
literature. Mathematically, MOG is complicated than GR in 
the sense that it postulates more gravitational filed than GR. 
In fact, MOG is a Scalar-Tensor-Vector theory of gravity, 
while GR is a tensor theory. In other words, in MOG in ad¬ 
dition to the metric tensor, there are two scalar fields (11 (x a ) 
and G{x a )) and also a massive Proca vector field (j) n . In fact 
the gravitational constant is variable as is in scalar-tensor 
theories of gravity. Therefore, the existence of these fields 
will provide some degrees of freedom which may help to 
handle the dark matter problem without invoking dark mat¬ 
ter particles. It is worth mentioning that adding non-minimal 
scalar fields to the effective Lagrangian of gravitational filed 
dates back to Brans-Dicke theory [2]. In this theory the non- 
minimal scalar filed was introduced in order to incorporate 
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the Mach’s principle into GR. However, there are also other 
motivations for adding into play the scalar fields. For exam¬ 
ple, non-minimally coupled terms between scalar fields and 
geometric quantities, such as the Ricci scalar R, have to be 
added to the gravitational action when quantum corrections 
are taken into account. Also scalar-tensor theories show up 
as the dimensionally reduced effective theories of higher di¬ 
mensional theories, such as Kaluza-Klein theory. For a com¬ 
prehensive review of the subject, we refer the reader to [3], 

On the other hand, the existence of the vector fields in 
the gravitational action has its own long story. In vector- 
tensor theories, in the form of Einstein-ather theories [4], 
there is a Lorentz-violating vector field. The existence of this 
vector field can dramatically affect the cosmology. For ex¬ 
ample, it may leave an imprint on perturbations in the early 
universe [5], and it can even affect the growth rate of struc¬ 
ture in the universe. Also the vector fields are a simple and 
natural candidates to explain some certain anomalies in the 
Cosmic Microwave Background [6]. For another example of 
theories containing vector fields, one may refer to Tensor- 
Vector-Scalar theory (TeVeS) [7]. TeVeS is a relativistic co¬ 
variant theory for Modified Newtonian Dynamics (MOND) 
[8], TeVeS has both types of above mentioned fields, i.e. 
scalar and vector fields. 

As we mentioned before, MOG is also another example 
of theories where scalar and vector fields participate in the 
gravitational sector. This theory has been applied to explain 
the rotation curves of spiral galaxies and the mass discrep¬ 
ancy in the galaxy clusters [9]-[ 14] . Recently, it has been 
claimed that this theory is in excellent agreement with the 
rotation curve data for the Milky Way, while MOND does 
not fit the data [15]. In fact, the observed rotation curve data 
for the Milky Way extends as far as 200 kpc from the center 
of the galaxy. Data for such a large distances from the core 
of the galaxy may provide a critical test of modified theories 
of gravity [15]. 
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There are some papers in the relevant literature which 
have studied the consequences of MOG at the cosmological 
scale, for example see [16]-[20]. However, to the best of our 
knowledge, the astrophysical aspects of this theory has been 
investigated more extensively than its cosmological conse¬ 
quences. For example, there is no exact cosmological solu¬ 
tion for this theory in the literature. This point motivated us 
to study the cosmological behavior of this theory by looking 
for exact solutions for the modified Friedmann equations. 
To do so, we use the Noether symmetry approach [21], In 
fact, as every conservation theorem, the Noether conserva¬ 
tion theorem allows to reduce dynamics of the cosmologi¬ 
cal model and gives insight into conserved quantities. This 
approach has been extensively used in various cosmologi¬ 
cal models and modified theories of gravity, for example see 
[22]-[33] . In [33] the Noether symmetry approach is used in 
non-minimal scalar-tensor theories and in higher order theo¬ 
ries of gravity. Also a general review of the Noether symme¬ 
try approach can be found in this paper. In the current paper, 
i.e. in the context of MOG, there is also a non-minimal scalar 
field. We shall see that there are some cosmological exact 
solutions in MOG which correspond to some specific solu¬ 
tions presented in [33], for example see subsection 6.4. It is 
worth mentioning at this point that the Hojman’s conserva¬ 
tion theorem [34] has also been used to study the dynamics 
of cosmological models and to find exact cosmological so¬ 
lutions for them [35]. 

The layout of the paper is the following. In Section (2), 
we briefly review MOG’s field equations and derive the cor¬ 
responding modified Friedmann equations. In section (4), 
we derive the point-like Lagrangian which is needed for ap¬ 
plying the Noether symmetry approach. Section (5) is de¬ 
voted to the application of the Noether conservation theo¬ 
rem to MOG. In this section we find the explicit form of the 
generic vector field potential by symmetry issues. In Sec¬ 
tion (6), for the first time in the literature, we find some new 
exact solutions for MOG and discuss their cosmological im¬ 
plications. Conclusions are drawn in Section (7). 


2 Modified gravity (MOG) 

We start with the generic action of MOG 

S=TLfV=gd 4 x[£(R-2A) + ^ v Vn X V V X d) 

A^-gP v S/ p y/N v V + ®o [\B^ v B^ v + V^]] + Sm 

where R is the Ricci scalar, Sm is the matter action, A is 
the cosmological constant, coo is a positive coupling con¬ 
stant, B pv = V^0 v — V v 0 /J , and the scalar fields % and l ft 
are related to those of the original paper [1] as % 2 =2/G 
and l// = Inp. Also is the self-interaction potential for 
the vector field 0“ and can be a function of 0 U 0 “ and iff. 
Variation of (1) with respect to g- u v , 0 a , l// leads to the 


following field equations respectively [36] 

Gyv +Ag pv = — j (V^ V v — g llv D)x~ + ~2 (2) 

A /C 


V^ + !!r = -— ja O) 

d<Pa OJo 

Ox=XR+^g tlv NvV'VvV' (4) 


□yr = 


“V r xvv+ 


2(0q dVfj, 

x 2 dv 


(5) 


Where G pv is the Einstein tensor and J a is a "fifth force" 
matter current defined as 


1 5S m 
\f~8 


( 6 ) 


nonzero J a means that there is a coupling between mat¬ 
ter and the Proca vector field 0“. Consequently, one can 
verify that the matter energy-momentum tensor is not con¬ 
served and the Einstein’s equivalence principle is violated 
[37]. Also, the total energy-momentum tensor is defined as 

nrtotal 7 - i _i _ (1\ 

1 yv - J fiv + + + \') 


where Tyv is the energy-momentum tensor for the ordinary 
matter, and 


T yv = (BfBva-gyv {\BP°B pa + V<>) + 2 ^) 
T yv = -jk ( V M* V vZ - \gyvVaXV a x) 

T yy = {VyYVvV'- 5 gyv^aV^ a v) 


Furthermore, we assume a perfect fluid energy-momentum 
tensor for the ordinary matter. Also we consider the cos¬ 
mological behavior of MOG in a flat Friedmann-Robertson- 
Walker (FRW) universe for which the line element is given 
by 

ds 2 = — dt 2 +a(t) 2 (dx 2 +dy 2 +dz 2 ) (8) 


where a(t) is the cosmic scale factor. Now, by inserting this 
metric to the field equations (2)-(5), and also noting the re¬ 
lation between j, y/ and G, /i, we find the modified version 
of the Friedmann equations, the 0-0 component of (2) gives 


8 tzG 



(9) 


G d 
Ga 


_1_ y 2 _i_ G 2 

12 j? 24 G? 


StcG 

3 


K(f)op 


Where Qo(t) = — <j>°(t) is the zeroth covariant component of 
the vector field and p is the total energy density of the mat¬ 
ter/radiation. Also dot denotes the derivative with respect to 
time. The terms inside the bracket can be considered as cor¬ 
rections to the corresponding standard Friedmann equation. 
Using the equation of motion of the test particles in MOG 
(the generalized version of the geodesic equation in met¬ 
ric theories), one can find an action for a pressure-less dust. 
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Therefore, it can be shown that J° = Kp (see [37] for more 
detail). Where K is another coupling constant. Although it 
seems that there are two independent coupling constants ©o 
and K, they almost always appear as the combination k 2 /coq 
through the calculations of the physical quantities. Also it 
should be mentioned that because of the cosmological prin¬ 
ciple, spatial components of J a and all components of anti¬ 
symmetric tensor B^ v are zero. Furthermore the i-i compo¬ 
nents of (2) takes the following form 


a 

a 


- 47 ^ } (p + 3p) + j+ ( 10 ) 

, r 1 ^ , 1 G 11 G 2 4tcG^a. n G«o v -i 

+ [ 2Ga + 6^ + 2 G-T2GZ-—^oP ~ ^} 


where p is the pressure of the cosmic fluid, and we have used 
equation (9) for terms including a 2 . Equation (3), i.e. the 
field equation of the vector field, takes the following form 


dVy, _ \6 kk 
<50o ®o 


( 11 ) 


Therefore the potential Vq should satisfy the following gen¬ 
eral criterion for all values of p(t) and 0o(f) 


<**> dv K > 0 
K <500 “ 


( 12 ) 


Furthermore, equation (4) can be written as 


G 

G 


327rG(l + K0o)p +12- +9 
a 



G d 
G a 
G 2 
G 2 


and finally equation (5) reads 


A 


p 2 pa 


Gp 
G p 


Gcoop 


dV± 

dp 


2G(tX)V(j, 


(13) 


(14) 


If the generic potential Vq and the equation of state p = 
p(p) are known, then equations (9), (10), (11), (13) and (14) 
are five non-linear partial differential equations for five un¬ 
known functions a(t), 0o(f), p(f), G(t) and p(t). It is inter¬ 
esting to remark that the vector field is not "dynamical" in 
the sense that its field equation does not contain any time 
derivative. In other words, the vector field is algebraically 
related to the energy density p. 


3 Vacuum exact solution 


Before starting the Noether symmetry approach, we con¬ 
sider the vacuum exact solution (p = 0 and A = 0). In this 
case the modified Friedman equations can be analytically 
integrated for a wide range of generic potential V^. Let us 
rewrite Equations (9), (10), (11), (13) and (14), respectively, 
as 


I _A 

12 p 2 + 2G 2 ’ 



(15) 


H 

77 


a_ 

G 


-3 


(16) 


« ft' 2 

p = -— 



(17) 

(18) 


where H = a/a and prime denotes derivative with respect to 
In a. Also it should be mentioned that we have assumed that 
the Proca potential has a general form V^> = f\ (p)/2(0 a 0“) 
where f\ and are arbitrary functions, satisfies / 2 ( 0 ) = 0 
and = 0. Note that in the absence of matter/radiation 

“TO 

the vector field freezes out at a finite value. We have as¬ 
sumed models for which the vector field freezes at 0o = 0. 
This is the case for original potential of MOG, i.e. Vq = 
-jP 2 0“0 a . 

The exact solutions to equations (15)-(18) are 


G(f) = Goa(t) e 


(19) 


±V 24e—24—£- 

p(t)=p 0 a{t ) v 72 (20) 

a(t) = (eg+ e 2 )3=s (21) 

where Go, po, £, £i and £2 are constants of integration. One 
may consider this solution as a late time solution where the 
energy density of matter and radiation is zero. If 2 < £ < 3 
then we have an accelerated expansion, i.e. a > 0. In this 
case G(t) is an increasing function of time and p (t) can be 
decreasing or increasing. It is important mentioning that this 
late time acceleration is not a de Sitter universe where the 
cosmic scale factor grows as a(t) ~ e'^ t . 


4 Canonical point-like Lagrangian for MOG 


Equations (9), (10), (11), (13) and (14) can also be deduced 
from a canonical point-like Lagrangian (a, a, %,%, 0o, 0 0 , VT V)- 
For notational simplicity, in this section we use y/ and % 
instead of p and G. The point-like Lagrangian .1 J is de¬ 
rived from action (1). We denote the configuration space by 
■2 = {a,X,VC0 0 }. Then = {a,d,x,X, 00,00, V, V} is 
the corresponding tangent space where the Lagrangian ££ is 
defined. On the other hand, the Euler-Lagrange equation of 
if is 

!§f-if = l (22) 

where q can be a, X or y/. Also the energy equation is 


<5if . dif . <5if . d5£ r/1 ,, 


(23) 


therefore combining equations (22) and (23), we have five 
equations which make a complete set of partial differential 
equations for five unknown functions. 
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In order to find Jz? one must insert the flat FRW metric 
into action ( 1 ) and simplify the result using an integration 
by parts. Note that the Ricci scalar curvature takes the form 
R = 6(a 2 /a 2 +a/a). As a result, one can easily verify that, 
in a FRW manifold, the Lagrangian related to the action (1) 
takes the point-like form 

Jz? = 6 Xiao 1 + 3ad 2 % 2 + Acc'i 1 + \a 2 ’x 2 + (24) 

^a 3 % 2 yf 2 - cooa^Vtj, + 16 npa 3 

where p = p m + p r is the proper energy density. Where p m 
is the matter energy density, and p, is the radiation energy 
density. It is worth mentioning that, in general, the mater 
energy-momentum tensor in MOG is not conserved (see equa¬ 
tion A31 in [37]). However in the case of FRW universe 
we have V U 7’’ UV =0 and consequently one can write p,„ = 
p,n0 a 3 and p,- = p^a 4 , where p,„o and p,-o are the current 
values of the corresponding energy densities. The current 
magnitude of the cosmic scale factor has been normalized 
to unity. It is interesting that although there is a coupling 
between matter and the vector field, the energy-momentum 
tensor of the ordinary matter is conserved. This is an enor¬ 
mous simplification made by the cosmological principle. 

Now by writing the Euler-Lagrange equations (22) for 
the Lagrangian (24), one can straightforwardly derive the 
equations (10)-(14). On the other hand, the energy condition 
(23) will give the equation (9). 


5 Noether symmetry approach 


Now let us introduce the lift vector field X [38]. X is an in¬ 
finitesimal generator of the Noether symmetry in the tangent 
space S7J2, and is defined as 


d d d d . d 

X = —f /3-j;—by^—b^ 3 — + 0 : 3 ^+ 

da d% dy/ <70o da 

& d . d , d 

+li Jx +r ^ + ^ 


(25) 


where the functions a, /), y and c depend on the configu¬ 
ration space variables a, X- V and 00 - One can straightfor¬ 
wardly show that 


— = L x Jz? + I6nt;a 3 Kp 
dt 


(26) 


where is the Lie derivative of the Lagrangian along the 
lift vector X, and E is defined as 


d££ <9 Jz? <9jz? 

E = a A^ + li W + r W 


(27) 


Therefor one may conclude that E is a conserved quantity 
if L x Jz? = — 16 tt %a 3 Kp (this is the Noether’s theorem). It is 
needed to recall that in theories which there is no coupling 
between matter and the gravitational fields, the right hand 


side of this criterion is zero. By applying this condition to 
the Lagrangian, the analytic form of the functions a, (3, y 
and 4 will be specified. In fact, one may simply equate to 
zero the coefficients of terms like y/ 2 , j 2 , y/f and so on in 
order to find some differential equations for the lift vector’s 
components. This is a common procedure in the relevant lit¬ 
erature of the Noether symmetry approach, for example see 
[21]-[33]. Applying this procedure to the Lagrangian (24), 
we find the following differential equations 

12zff + 12«§f+Z« 2 ^ = 0 

P a +Tl£ + X 2 jl +a X $+ax^+2a X = 0 

|a + 6 *ff+ag =0 (28) 

Xa+2ax^ +2a[i +2 a 2 ^ = 0 
3ax+2ax^+2ap=Q 
12 zff+ 2 a^+«Z 2 g =0 

dy _ da_ _ df_ _ df_ _ Q 
dtjo d tj> 0 d(j)Q dj> 0 

Lurthermore, there are some terms in which do not 

contain any time derivatives of configuration space variables. 
The summation of these terms should be separately zero: 


(3co Q V l j ) + 167ip r )a-2Axa + 

dV 

cooa^y= 0 

dy/ 


(29) 


It is obvious from equations (28) that the function ^ can be 
an arbitrary function of a , x an d y/. In other words, the ex¬ 
istence of the Noether symmetry does not restrict the func¬ 
tional form of 4 . Therefore, one can construct infinite num¬ 
ber of symmetry generators for this cosmological model. 
It is also worth mentioning that, there are nine differential 
equations in (28) which are enough for finding the exact 
form of ot(a,x, yt), y/), y(a,X, V / )- After finding this 

functions and inserting them into (29), we shall find the form 
of the generic potentials for which the Noether symme¬ 
try exists. After some straightforward algebraic calculations, 
one can verify that the general solutions to equation (28) are 


a = (aiy/ + a 2 )a 
y = 6a\ In (ax~ 3/2 ) + «3 


(30) 


Where ofi, cti and 0:3 are integration constants. Although 
the Lagrangian (24) seems complicated, it possesses simple 
Noether symmetry generators. Lurthermore, using the sec¬ 
ond equation in (30), the term containing the cosmological 
constant A vanishes in (29). Therefore, equation (29) can be 
rewritten as 


V(/) + 


\6np r 
3 <X>0 


6«i In (ax 3/2 ) + OS3 dV$ _ Q 
3(a\y/ + a, 2 ) dy/ 


(31) 
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In order to find the exact form of the potential Vg it is needed 
to find the scale factor a(t) and the scalar field % with respect 
to y/ and 0o- To do so, let us first derive the conserved quan¬ 
tity E associated to the Noether symmetry generator X. We 
recall that Noether’s theorem states that every symmetry of 
the Lagrangian implies the existence of a conservation law, 
and therefore a conserved quantity. Using the Lagrangian 
(24) and equation (27), we find 

E = (9a 2 xjC-6x 2 aa)a+a i X 2 ¥Y (32) 


For the sake of simplicity we assume that E = 0. Thus one 
can easily verify that when ci\ ^ 0 , we have 

3 OC 3 

ln(fl^ _ 2 ) = c^ai - - -b 0 : 40 : 1 ^ (33) 

6(X\ 

where CC 4 is a new integration constant. When GC\ is zero (i.e. 
when 7 is constant), it is easy to show that 

/ 3. cc A 

In {ax J) = a 4 + — yr (34) 

oai 

Substituting these results into equation (31), we find 


Vg + 


167rp r 
3 (Do 


+A - 1 



(35) 


where A is a constant and is equal to when oq 0 and 
A = ( 3 a 2 )/ct 3 when oq is zero. It is important to note that, in 
principle, one may find the cosmic scale factor as a function 

dVfi 

of from the field equation of the vector field (i.e. equa¬ 
tion (11)). Therefore, equation (35) is a partial differential 
equation for potential Vg(tj) o,V / )- This equation will deter¬ 
mine the general form of Vg for which the Noether symme¬ 
try exists. 


5.1 The case p r = 0 and p m , A ^ 0 

In this case the energy budget consists of the cosmological 
constant A, non-relativistic matter, the scalar fields G and p, 
and the vector field 0- u . Equation (35) takes a simple form 
and can be easily integrated. The result is 

Vg(p,<$> 0 )=f(fa)p~ A (36) 

where /( 0 o) is an arbitrary function of 00 , and we have 
used y/(t) = ln(ju(f)). This potential covers a wide range of 
potentials. For example the original potential of MOG, i.e. 
Vg = — \p 2 (p v (f> v [1], lies in this subclass. Another, straight¬ 
forward example is 

V^p,r) = V 0 p n (-M v )' n/2 (37) 

which in an isotropic and homogeneous space-time takes 
the form Vg = Vop n (j) Also, it is obvious that if ao= 0 (or 
equivalently A = 0) then Vg is a function of the vector field 
and does not depend on the scalar field p. More importantly, 
every function of 0 ^ 0 ^ is a solution. 


5.2 The case p m = 0 and p r , A ^ 1 


In this case, we rewrite equation (35) as follows 
Vg 


' gt +A -'*i = 0 

3k dtpo dy/ 


(38) 


and the general form of the solution is 

Vg=f(y)p- A (39) 

in which y = exp( 0 o)p~ a / 3k: , and / is an arbitrary function 
of y. Simple examples of this class of potentials, are Vg = 
V 0 p and Vg = V o exp(A lx /-0 a 0“)p ;i2 . 


6 Cosmological exact solutions for MOG 

Now we present some new exact cosmological solutions to 
MOG’s modified Friedmann equations. Also we restrict our¬ 
selves to solutions with zero radiation density. We recall that 
by setting the Noether constant of motion E to zero, we 
found equations (33) and (34). Let us rewrite these equa¬ 
tions as 

a(t) = cG(f) - 3 / 4 ;U(f) 1/2A (40) 

where c is a constant. This is a key equation to derive ex¬ 
act solutions. Using this equation we find some solutions for 
specific potentials introduced in the previous sections as po¬ 
tentials for which the Noether symmetry exists. 


6.1 Model Vg = Xop{t) 2q , p m = p r = A = 0 


Where Ao is a constant. We have used a new parameter q 
instead of A as A = —2 q. This model corresponds to a scalar- 
tensor theory with two scalar fields p{t) and G(t). In fact, 
since the matter and radiation are absent in the model, there 
is no source for producing the vector field. With the aid of 
equation (40) and assuming a power-law type solutions for 
the scalar fields as p{t) = pot" and G{t) = Got"', we find the 
following exact solution 
1 

a(t) =a 0 t *- 2O ‘< 2 

G(t) = G0t^A (41 ) 

p(t) = Pot'-Si 1 


where Go is related to other constants as 

5 20 q 2 +7 

8 (Oq Ao po 2q (25 g 4 — 10g 2 + 1) 


(42) 


One can easily show that this solution corresponds to an ac¬ 
celerated expansion (i.e. a > 0 and d > 0 ) if \J 3/20 < q < 
\/y/~5 or — 1 /x/5 < q < —^3/20. In this case, the mag¬ 
nitude of the scalar G(t) decreases with time and for p(t ) 
increases. 
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This solution has an interesting feature: the model can 
be contracting (i.e. a < 0) while a > 0. This situation hap¬ 
pens when q > 1/V5 or q < — \/y/5. In this case p(f) is a 
decreasing function of time and G(f) is increasing. In fact 
the existence of the scalar fields p(t) and G(t) can lead to 
repulsive gravitational force in the early universe. 


6.2 Model = — ^p 2 (j) a (l) a , p r = A= 0 and p m ^ 0 


As we have already mentioned, this potential is the origi¬ 
nal potential postulated in MOG [1]. For this model we use 
equation (40) and also assume power-law type solution for 
G(t) and p(t ). One can easily verify the following exact so¬ 
lution 



544 71 CCoGfi/PmO , 
9 Mo 2 r 


12 

23 


G(t ) 


_ 2312 oiqGn 

- 207 ^2 


/X - 18 

p(t) = lM)t H 


10 

2 J 


= It 


Pm (t ) 


= _9_ Mg - 
544 tcocqGm 


36 

23 


(43) 


where olq=k 2 /(OqGn is one of the free parameters of MOG, 
and G n is the gravitational constant. Its observational values 
is gcq = 8.89 ± 0.34 [12], One can easily verify that the con¬ 
dition (12) is satisfied. It is obvious that the "velocity" of 
the vector field is zero and field’s value remains constant in 
this matter dominated phase. Also this solution corresponds 
a decelerating universe. If we compare this solution with the 
matter dominated universe in GR (flat FRW space time with 
A = 0) where a(t) ~ f 2 / 3 , then we see that the scale factor 
increases slower than the standard Friedman model. Thus 
the matter density also decreases with slower rate. 


6.3 Model = Vo p Pr = A =0 and p m ^ 0 


Where Vo is a constant. With the aid of equation (40), one 
may easily verify the following solution 


a(t) = aot 
G(t) = G 0 r 2 ~ m 

,.( t \ _ K 2 (3nm+10 k) 2 n 
' 4Go^ 2 V 0 hi(Bo 

K Mt) = - (3^,+lOzr) 
„ _ 3 jrm + 10 jr,m 

- 16G 0 7i- 1 


(44) 


where m can has two values: m = —18/7 and m = —3, and 
Go is related to the other constants as 


G 0 = 


(3nm + 107r)oo 3 

16 PmO ft 2 


(45) 


it is worth mentioning that the condition (12) holds if Go > 
0. As we see the vector field is again constant. Also, the exis¬ 
tence of two different values for m shows that our power-law 
type solutions are not the complete and unique solutions to 
the cosmological field equations. For m = —3 the cosmic 
scale factor increases linearly with time and so without de¬ 
celeration, i.e. a = 0. therefore scale factor increases faster 
than the standard matter dominated Friedmann model. Also 
we recall that in GR, for a matter dominated universe with 
flat spatial curvature and zero cosmological constant, a can 
not vanishes. 

For m = —18/7 the scale factor increases as a{t ) ~ f 6 / 7 . 
therefore the scale factor increases faster than the standard 
case. However, unlike the case m = —3, we have a negative 
acceleration (a < 0). As a final remark for this solution, we 
mention that the repulsive nature of the theory is again ev¬ 
ident. In other words, the combination of the scalar fields 
and the vector field MOG can yield a repulsive gravitational 
force at the cosmological scale. 


6.4 Model = Vo v 7 — (j> a (l> a p'\ p, = A = 0 and p m ^ 0 


Where m and n are real constant, and m / 0,1,2. Therefore 
this model is different from those presented in the previous 
subsections. Again, by using equation (40) which was ob¬ 
tained from the conserved quantity of the Noether symme¬ 
try, we find the following solution 


a[t) = aot 2 

G(t) = Got 4 

p(t)=p 0 t K 

K< h(t) = K&0 = JpAmj 

O ( f ) — 5 ( 2 — m ) 
PmVJ— 8 jrGo t 6 


(46) 


where K = ±2y/l5 corresponding to n = ^ \/l 5/5 and the 
relations between constants are 

10p ( f /IO <V" 5(2 -m)al 

Vo = r„r ’ Go= < 4? ) 

(OoGo on p m o 

The criterion (12) leads to (2 — m)/Go > 0. If we assume 
that Go >0 (note that G(t) is related to the gravitational con¬ 
stant), then everything is well defined provided that m < 2. It 
is interesting that this solution represent a constant acceler¬ 
ation for the cosmic expansion. We mention again that such 
a behavior is not achievable in the corresponding standard 
Friedmann model. As it is clear from equations (46), the 
scalar field G(t) is an increasing function of time, and the 
scalar field p(t) depending on the value of parameter K can 
be an increasing function of time or a decreasing one. How¬ 
ever, the rate of expansion is independent of the behavior of 
P(t). 


>0 
















7 


It is worth mentioning that G is a non-minimal scalar 
field. As we have already mentioned, non-minimal scalar- 
tensor theories have been already investigated using the Noether 
symmetry approach, for example see [33], Although there is 
an extra vector field in MOG which is coupled to matter 
distribution, there are some similarities between the exact 
solutions of these theories. For example the exact solution 
presented in this subsection is reminiscent of the exact solu¬ 
tion presented in equation (6.42) in [33]. 

6.5 Model Vt = ~\n 2 <l> a <t>a + ^V^F r ^+£, 
p r =A= 0 and p m ^ 0 

In this case, we found the solution exact solution: 
a(t) =a 0 f- 3 D m+ i 
G(t) = Got t 

p(t) = p 0 t' n (48) 

k:< j>o(t) = -1 - 

V 10 Go Mo 2 

a o 

where m = 5/4 or —20/II and 

ao = ?V23 2 /3 5 V6 s/ GoaoG^W 
3 y m z p o 

The condition (12) is satisfied if po > 0 and Go > 0. For 
m = 5/4 the scale factor decreases as a(t) ~ /~ 9 / 8 , and so 
a < 0 and a > 0. However, for m = —20/11 the scale factor 
grows as a(t) ~ t. 

7 Conclusion 

In this paper we have studied the cosmology of a specific 
scalar-vector-tensor theory of gravity known as MOG using 
a well-known approach. We have found the Noether sym¬ 
metry generators and also the conserved quantity associated 
to the symmetry. More specifically, this approach helped 
us to find the explicit form of the vector field potential Vy> 
(see equations (36) and (39)). We showed that the origi¬ 
nal potential postulated in this theory, is one of those po¬ 
tentials for which the Noether symmetry exists. Also, using 
the conserved quantity associated to the Noether symmetry, 
we have found some exact solutions for the modified Fried¬ 
mann equations of MOG. We have found five set of exact 
solutions and except in one of them, we have not neglected 
the ordinary matter contribution. More specifically, our solu¬ 
tions correspond to matter dominated universes. As we men¬ 
tioned before, power-law type exact solutions are not neces¬ 
sarily unique. However, these solutions are very often used 
in cosmological theories in order to insight into the physical 


content of them. In the context of MOG, there are some in¬ 
teresting features in the exact solutions which may help to 
understand the cosmological behavior of this theory. For ex¬ 
ample, the existence of the scalar and vector fields in MOG 
can lead to contracting universe with positive acceleration. 
In other words, these fields effect can behave like repulsive 
gravitational force at the cosmological scales. Also, there 
exist some solutions which the "gravitational constant" (the 
scalar field G(t)) increases with time as well as solutions 
where G(t) is decreasing. 
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